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Abstract 



We show the way in which the self-consistent Ornstein-Zernike approach 
(SCOZA) to obtaining structure factors and thermodynamics for Hamil- 
tonian models can best be applied to two-dimensional systems such as 
monolayer films. We use the nearest-neighbor lattice gas on a square 
lattice as an illustrative example. 



I. INTRODUCTION 



The self-consistent Ornstein-Zernike approach (SCOZA) was introduced by H0ye and 
Stell as an approximation method specifically tailored to obtain structure factors and 
thermodynamics for Hamiltonian models in three or more spatial dimensions. It was sub- 
sequently found by H0ye and Borge |2| that the SCOZA yields extremely accurate results 
for the two-dimensional lattice gas as well, when appropriately used, thus opening the way 
toward the use of SCOZA in treating thin-film problems. In this article we summarize the 
two-dimensional SCOZA results. We point out that those results for systems considered 
in the thermodynamic limit are strikingly similar to the results that would be found in an 
exact analysis of two-dimensional systems that are finite || or semi-infinite 0. We note 
why this is to be expected, and using the behavior of the specific heat as a criterion, we 
find the size of the finite and semi-infinite systems that yield the best match to the SCOZA 
results for the infinite square lattice. 



II. BACKGROUND 

The SCOZA is based on an ansatz used by Ornstein and Zernike ||, which is that 
the direct correlation function c(r) introduced by those authors has the range of the pair 
potential. In SCOZA, this ansatz is used along with a core condition that guarantees that 
the two-body distribution function g(r) must be zero for values of r for which the pair 
potential is infinite. In a lattice gas, the core-condition implies the exclusion of multiple 
occupancy of a single lattice site or cell; for the equivalent Ising model in which the spin 
variable at site i is +1 or -1, the corresponding condition is that (siSj)i=j = 1, which simply 
reflects the fact that the spin must be pointing either up or down with probability one. 

An analysis of the Ornstein-Zernike formalism made by one of the authors some time 
ago showed that the core condition plus the assumption that c(r) is proportional to the 
pair potential implies that for short-ranged potentials, a two-dimensional system can not 
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have a critical point at nonzero temperature. (In three or more dimensions there is no 
such restriction on criticality.) One knows however, that in two dimensions, systems such 
as the nearest-neighbor lattice gas do in fact have a critical point at nonzero T c . Hence 
SCOZA did not initially appear to be a promising method for treating arbitrarily large two- 
dimensionall systems. Its apparent unsuitability is also consistent with the observation |||7j 
that the assumption that c(r) has the range of the potential implies that for short-ranged 
potentials, the critical exponent 77 = 0, so that at a critical point at T c 7^ 0, where on expects 
g(r) — 1 « const ./ r d ~ 2+v , SCOZA would yield g(r) — 1 ~ const. /r d ~ 2 . In three dimensions, 
in which n « 0.03, and d — 2 + n « 1.03, this leads to negligible error. But for d = 2, in 
which 77 = 1/4, it represents the difference between g(r) — 1 ~ const. /r 1 / 4 and a g(r) — 1 
that does not appropriately decay with increasing r. 

However, as we shall see below, in two dimensions, the SCOZA results for a square lattice 
of infinite extent are strikingly similar in some respects to exact results for either an N x N 
lattice or an N x 00 lattice, with N « 22. 

It is not hard to understand why the assumption of a c(r) for an infinite square lattice 
yields results that mimic exact results for a finite system. For a finite system, c(r) is limited 
in range by the finite boundaries of the system. One also knows that in an exact analysis of 
an infinite one- dimensional system with short-ranged potential one finds no critical behavior 
for non-zero temperature. From these qualitative statements however, it is not clear what 
values of N in an iV x N or iV x 00 lattice will give rise to exact results that most closely 
match SCOZA results for an infinite square lattice. As we shall see in Section IV, when one 
uses the behavior of the specific heat as a criterion, N turns out to be around 22. 

III. THEORY 

In the following we consider the two-dimensional square lattice gas, which is isomorphic 
to the two-dimensional Ising model. The potential between particles is 
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v(ri — rj) — <i — -u;^ i ; j nearest neighbors (1) 
0, otherwise. 

For convenience w is scaled to be 1 in our calculations. In this convention the internal energy 
per spin for the Ising model U is related to the internal energy per particle for the lattice 
gas u via the following relation: 

1 1 

U = pu + -qp--q, (2) 

where q is the number of nearest neighbors (q — 4 for the square lattice), and p is the density 
for the lattice gas. 

SCOZA is based on the enforcement of thermodynamic consistency between different 
routes to thermodynamics. This imposes the following relation: 

d(3 dp 2 ' 1 ' 

where f3 = 1/T, the inverse temperature, p~ 2 x is the isothermal compressibility, and u is 
the internal energy per particle for the lattice gas. These quantities can be given in terms 
of correlation functions through fluctuation theory, which yields 

PX' 1 = — , (4) 

p(i + pM0)) 

and through the ensemble average of the Hamiltonian, which yields 

« = -^QP9i = ~QP0- + hi) , (5) 

where h(r) = g(r) — 1, and g(r) is the two-body distribution function. Here g± and hi 
represent the functional values of g(r) and h(r) at nearest-neighbor positions, and h(k) is 
the Fourier transform of h(r), which is related to the direct correlation function c(r) by the 
Ornstein-Zernike equation: 

h(r l )=c(r l )+p'£c(r J )h(r l -T J ), (6) 



or in the Fourier-transformed space: 

1+ "^ = r^(k)- m 

The above relations are exact. In order to proceed we shall approximate the form of the 
direct correlation function c(r) by using the ansatz introduced by Ornstein and Zernike (OZ) 
that c(r) has the range of the pair potential. In SCOZA this can be done by generalizing 
somewhat the mean spherical approximation (MSA), in which 

c(k) = c + gci$(k) , (8) 

where $(k) is the nearest neighbor sum, 

,_ v cos kx ~ \~ cos ky . , . 

<P(k) = , tor a square lattice, (9) 



and c and Ci are functions of (p,/3). Eq. (g) is the OZ ansatz applied to the lattice gas. 
The MSA is the special case obtained by setting C\ — f3 and adjusting Co to be compatible 
with the core condition that assigns zero probability to multiple occupancy of a single site. 
In SCOZA one instead adjusts C\ to insure self consistency between Eq. (f|) and (H). 

The core condition h(0) = — 1 implies a relation between Co and c\ through the OZ 
equation: 

r d 2 k 1 
1 ~ P ~ J (2^)2 1 - pc(k) 

1 r d 2 k 1 



1 - pc J (2vr) 2 1 - z$(k) 

./Si, do) 

1 - pc 

where z = qpc\/ (1 — pc ). P (z) is the value for the lattice Green function P(z, r) at r = 0. 
For a two-dimensional square lattice we have 

P { z) = 1 / f - ± = ljf W , (11) 

where i^(z) is the complete elliptic integral of the first kind. From Eq. (|10D we have 



Co 



p 



p{zy 

1-p 



and 



Cl 



gp(l - p) 



By taking r ; = in Eq. @, we get 



- 1 = /i(0) = c - pc + gpci/ii 



(12) 



(13) 



(14) 



so 



1 . , N . 1 - p 1 - P(z) 

1 + 1 - p c = g — y 

gpci p zP(z) 

After substitutions Eq. (§) and @ become 

:i-^)pu) 



p(l - p) 



and 



1 ( n n 1 "^) 



(15) 



(16) 



(17) 



From Eq. ( jl6|) we conclude that in SCOZA the criticality, if any, occurs at z = 1, since at 
the critical point one has x" 1 — 0, and we have P(z) > for all z. 

Finally, by applying the thermodynamic consistency via Eq. (|3p, we get the SCOZA 
partial differential equation: 







p(l - p) d(3 



[(l-z)P(z)] = - 



q d 2 
28^ 



P(l-P) 



zP(z) 



-q- 



(18) 



The boundary conditions are z = 0, i.e., P(z) — 1, at (3 = and p = 0, 1. 



IV. SPECIFIC RESULTS AND DISCUSSION 



Since for the two-dimensional square lattice P(z) diverges when z — > 1, the renormalized 
inverse temperature parameter ci also diverges. As a result, SCOZA fails to predict a true 
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critical point above zero temperature in this case. Nevertheless, the SCOZA results for u 
match the exact Onsager expression || for u remarkably well over the whole temperature 
range. Instead of having an infinite slope at the exact critical temperature, the SCOZA slope 
achieves its maximum at a temperature within a fraction of a percent of the exact value. 
The nonsingular but near-singular behavior near the ideal transition temperature makes our 
results strikingly similar to the exact results for a finite-size Ising model on a square lattice, 
iV x iV ||, or a finite- width strip, iV x oo |], for an N a bit greater than 20. 

In Fig. 1 we plot the negative internal energy —U versus the inverse temperature (3 along 
the critical isochore p = p c = 1/2, i.e., magnetization being equal to 0, for the comparison 
between SCOZA and both infinite- and finite-size Ising exact results. The comparison 
is made even clearer by plotting their residuals in Fig. 2. We find that the deviations 
between SCOZA and the other results are very small over the whole temperature range. 
The deviations get larger near the critical point (3 = /3 C , although the largest deviation is 
still within 3 percent. Comparing with finite-size and finite-width exact results, we find this 
deviation gets minimized when we choose N = 22 for an Af x iV Ising model or N = 21 for 
an N x oo Ising model. 

In Fig. 3 we plot the specific heat versus (3 along the critical isochore. The SCOZA result 
has a specific heat that stays finite at its maximum, as all the finite-size and finite-width 
specific heats do in an exact theory. In this comparison we again find that the SCOZA 
infinite-lattice result has great resemblance to the 22 x 22 or 21 x oo Ising model. The 
SCOZA curve matches well with the 21 x oo Ising model for (3 < (3 C . But for (3 > (3 C 
the SCOZA result has somewhat less difference to the 22 x 22 solution. The maximum of 
SCOZA curve occurs at (3scoza = 1-758, whereas the exact critical point for the infinite 
Ising model is (3 C = 1.763. The maximum for the 22 x 22 and 21 x oo Ising models occur 
at /?22x22 — 1-735 and /32ixoo = 1.764, respectively. The deviations between these maximum 
temperatures and the exact critical temperature is of the same order for the SCOZA and 
21 x oo results, whereas it is a bit larger for the 22 x 22 case. In this regard the SCOZA is 
a bit more similar to the 21 x oo model. 
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For each fixed value of p, we define the temperature where the specific heat is at its 
maximum as the pseudo singular temperature. Here we try to determine a 'pseudo' spinodal 
curve by collecting the set of (p, T) points that correspond to those pseudo singularities. This 
'pseudo' spinodal curve is shown in Fig. 4. Note that the points on the curve don't mark 
real singularities, and both the specific heat and the compressibility remain finite, but very 
large, at these points. Furthermore, near the pseudo critical point p = p c = 1/2, by defining 
Ap = p — p c , AT = T c — T, we find 

Ap ~ (AT)^" 1 ^ , (19) 

with a classical exponent (5 sp i no dai = 1/2 when T < T c . 

From Eqs. (|7|) and (^j) we find that in SCOZA the function 1 + ph(k) has a form pro- 
portional to the Fourier transform of the lattice Green function P(z, k). Hence we have 
S(r) + ph(r) oc P(z,r). From the fact that 

h(v) ~ P(z, r) ~ exp(-2r^^) (20) 

when r — > oo |§, we find in SCOZA the correlation length 

H/S- < 21 > 

In Fig. 5 we plot the correlation length £ versus temperature at the critical isochore. Since for 
the two-dimensional SCOZA scheme z — ► 1 only when T — > 0, the correlation length keeps 
finite at the exact critical temperature T c and only diverges when T — ► 0. However, when 
T = T c we have £ ~ 44 according to the SCOZA scheme, which already indicates strong 
correlations. Furthermore, the correlation length increases sharply right below T = T c . 
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FIGURES 

FIG. 1. Negative internal energy for the two-dimensional Ising model —U, versus the inverse 
temperature (3. The solid curve is the SCOZA result. Other curves are exact solutions for 
22 x 22,21 x oo, and oo x oo Ising models, respectively. (3 C is the exact critical point for the 
co x oo model. 

FIG. 2. Residuals between SCOZA and exact results. For optimal finite-size results, the devi- 
ations are smaller than the oo x oo exact results. 

FIG. 3. Constant volume specific heat Cy for SCOZA compared with exact results. Note that 
the specific heat for SCOZA doesn't diverge at its maximum, showing a resemblance with finite-size 
models. 

FIG. 4. Pseudo spinodal curve in the (T, p) plane, derived from SCOZA. Near the pseudo 
critical point p = 1/2, there is the relation Ap ~ (AT)^™ * 2 ' with /3 sp i no dai = 1/2. 

FIG. 5. Correlation length £ versus temperature from the SCOZA results. Note that £ starts 
to increase sharply for T ~ T c . 
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